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I. INTRODUCTION 

Recent progress in Hall experiment with graphene reveals many new features in longitudinal and transversal re¬ 
sistivity in Hall configurations characteristic for fractional quantum Hall effect (FQHE) both in suspended graphene 
scrapings EHi, and in graphene samples on crystalline substrate of boron nitride [3, 0 • The different structure of 
Landau levels (LLs) in graphene in comparison to the conventional semiconductor 2DEG is the source of distinct 
scheme for integer quantum Hall effect (IQHE) in graphene referred as to ’relativistic’ its version Q. The Berry 
phase induced shift for chiral carriers in graphene together with the four-fold spin-valley degeneration of LLs result in 
;/ = 4 (n -b i) series for fillings at which IQHE plateaus occur in subsequent centers of LLs for the case of monolayer 
graphene. When SU{A) spin-valley symmetry is broken by larger magnetic field, the new features for IQHE appear 
corresponding to removing of LL subbands degeneration Q. In bilayer graphene the extra degeneration of n = 0 
and n = 1 LLs shifts plateau positions to ends of LLs also four-fold degenerated in this case except for the eight-fold 
degenerated LLL. Simultaneously, still more and more features at fractional fillings of LLs are observed related to 
FQHE, revealing also a specific its character mostly connected with the different structure of subbands of LLs and 
particle-hole symmetry due to meeting of the conductivity and valence bands in Dirac points. The new filling frac¬ 
tions are observed in six first subbands of LLs with n = 0 and n = 1 in monolayer graphene not repeating the 
hierarchy of FQHE in conventional 2DEG and not allowing for explanation upon standard composite fermion model. 
Especially interesting is an observation of unusual even denominator Filings for FQHE in bilayer graphene including 
the most pronounced feature at iz = — i [3|, taking into account possibility for control of non-Abelian anyons excited 
from the — ^ state as previously studied for = | in traditional semiconductor 2DEG Q and important for potential 
application in topological quantum information processing. 

In the present letter we analyze the details of the hierarchy for fractional Filings linked to strongly correlated 
multiparticle states in graphene using the topological commensurability approach developed earlier for the ordinary 
2DEG Hall systems |9l-Hll| . In this way we explain the structure of fractional Filings of LL subbands and demonstrate 
its evolution with growing number of the LL and the number of its subband. This approach gives the hierarchy of 
FQHE and of other features, including paired states and Hall metal, in agreement with the available experimental data 
for monolayer and bilayer graphene. In particular the explanation of the even denominator Filing ratios for bilayer 
graphene, with the FQHE state at = — i, has been achieved in a direct manner upon the applied commensurability 
topological method. 

II. COMMENSURABILITY CONDITION—THE BASIC FORMULATION 

The concept of commensurability of classical cyclotron trajectories with interparticle spacing in 2D systems is born 
in relation to the braid group approach to multiparticle systems in the presence of a magnetic Feld. If N identical 
particles are located on the manifold M, the collective behavior of such system can be assigned by the statistics 
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phase shift if one considers position exchange of particle pair I, IS- This quantum feature is associated with the 
one dimensional unitary representation (IDUR) of the full braid group related to the system. The full braid group 
is the first homotopy group of the multiparticle configuration space, 7ri(<l>) = — A)/Sn), where is the 

iV-fold normal product of the manifold M, A is the diagonal point set in this product (when coordinates of two or 
more particles coincide, subtracted in order to assure conservation of the number of particles), S'tv is the permutation 
group of N elements [l^- The quotient structure of the configuration space $ = {M^ — A)/Sn) is related with the 
indistinguishability of quantum identical particles—the property essential for the quantum statistics determination 
reason. Apart of this quantum prerequisite the full braid group 7ri($) is the classical topological object collecting all 
classes of nonhomotopic trajectory loops in the configuration space $, where points which differ only by enumeration 
of particles are unified. Any details of the dynamics of N particle system caused by the interaction and resulting in 
special shapes of trajectories are not important here. Only topology of trajectories decides whether one trajectory 
loop can be continuously transformed into another one or not. In the latter case such topologically nonequivalent 
trajectories fall to the distinct classes of the full braid group. Thus, the full braid group does not reflect the dynamics 
details but rather identifies the topology restrictions imposed on the multiparticle system associated to some global 
features of the system including the manifold M type. As the loops from 7ri($) describe exchanges of particles 
(due to identification of positions of all particles which differ by the enumeration only) the full braid group contains 
information on quantum statistics of particles, though classical particles do not have any such statistics. It was 
proved that the one dimensional unitary representations (IDURs) of the full braid group serve as determinants of 
the quantum statistics of quantum particles related to the original classical ones [l^. Therefore, if particle classical 
positions defined by the coordinates of the multiparticle wave function 'I'(xi,..., xvr) {xi is the coordinate of i-th 
particle on the manifold M, i.e., the classical position of the i-th particle on M) are changing along a selected loop 
from the 7ri($), then this wave function acquires the phase shift e®“ given by the IDUR of this particular braid [l^. In 
this way the statistics of quantum particles can be identified. For the same classical particles the quantumly different 
particles can be defined as assigned by distinct IDURs of the related full braid group. For three dimensional (or 
of higher dimension) manifolds M the full braid group is simply the permutation group S'at, with only two different 
IDURs [13, 


( 1 ) 

where Gi is the generator of Sn, i-e., the elementary braid describing exchange of i-th and i -I- 1-th particles, whereas 
other particles are left in their positions. These two IDURs correspond to bosons and fermions, respectively. For two- 
dimensional manifolds M the braid groups differ considerably from the permutation group. Unlike the permutation 
group, for M = R^, 7ri($) is an infinite group with the IDURs (l3 - 

CT* as [0,27r), (2) 

corresponding to so-called anyons (including bosons for a = 0 and fermions for a = tt). Anyons are quantum particles, 
different than bosons or fermions, which can appear in 2D multiparticle systems, e.g., for particles located on the R? 
plane or on the locally 2D manifold, like a sphere (but not on torus [l3, for which IDURs of related full braid group 
do not exist). 

Nevertheless, there are also other peculiarities of 2D manifold topology, which were not accounted for in the 
manner as described above. The special feature of planar multiparticle systems manifests itself in the presence 
of the perpendicular magnetic field strong enough that the classical cyclotron radius is shorter in comparison to 
the interparticle separation on the plane. Because the classical trajectories of free charged particles are defined 
by cyclotron orbits at the presence of a magnetic field, thus exchanges of neighboring particles in 2D topology are 
possible only if the size of the cyclotron orbits fits to the separation between particles, as is illustrated in Fig. [T] Let 
us emphasize that the average separation between particles on the plane is rigidly fixed by the Coulomb interaction 
between particles preventing approaching one particle onto another one in the uniformly distributed N particles. The 
fulfillment of the commensurability between the cyclotron radius and the interparticle spacing is required for the 
definition of the generators Ui of the full braid group, where Gi describes exchange of i-th and i -I- 1-th particles. If the 
cyclotron orbit is incommensurate with the separation between particles on the plane, then the classical trajectories 
describing Gi elements of the full braid group are impossible and this braid group cannot be defined for such strong 
magnetic field. It means that the quantum statistics cannot be defined in this case and any correlated multiparticle 
quantum state cannot be organized. This situation happens at fractional filling of the lowest Landau level (LLL). 
For magnetic fields larger that this one which corresponds to the completely filled LLL, the classical cyclotron orbits 
are too short to match neighboring particles and the definition of the generators of the full braid group is precluded. 
Nevertheless, at some magic fractional filling ratios of the LLL the correlated states are experimentally observed and 
refereed as to fractional quantum Hall effect (l6l-[l8l|- In means that possibility of particle exchanges is recovered 
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somehow. In the framework of the composite fermion (CF) approach [l^ the enlargement of cyclotron orbit size is 
achieved by screening of the external magnetic field by auxiliary field flux quanta attached to particles. This artificial 
model allows for identification of the main line of fractional fillings for FQHE by mapping of the fractional state onto 
integer quantum Hall states in resultant magnetic field reduced by the average field of fluxes pinned to CFs. This 
model does not explain neither the origin of auxiliary field fluxes nor the mechanism of creation of composite particles. 

The manifestation of FQHE can be explain, however, also in the braid group terms 0- Though the generators at of 
the full braid group cannot be defined when the cyclotron orbits are shorter than interparticle separation, there exist 
other braids which in 2D fit to interparticle separations. This exceptional property of planar multiparticle systems 
possess multilooped cyclotron braids, i.e., where q-odd integer, and for e.g., q = 3, af describes the braid for 
exchange of z-th and i + 1-th particles with an additional loop. Exclusively in 2D the multilooped cyclotron orbits 
have larger size which can fit to interparticle separation at the magic filling fractions of the LLL, the same ones at 
which FQHE is observed. The reason of the enhancement of the size of planar multilooped orbits is linked with 
constant surface field of 2D orbits despite its multilooped character (in opposition to 3D case when each additional 
loop adds also a surface portion spanned by this loop). When the total external field flux is passing through the 2D 
multilooped orbit it must be shared between all loops and their sizes grow. This is illustrated in Fig. [2] 

In Fig. [2] (left) the scheme of the cyclotron orbit at magnetic field B is shown as accommodated to the quantum 
of the magnetic field flux, i.e., BA = This serves here as the definition of the cyclotron orbit A because this A 
fits to the interparticle separation, S —the sample area, N —the number of particles, in the case of the completely 

filled LLL. If only single-loop orbits are available, then at, for instance, 3-times larger field, 35, the cyclotron orbit 
accommodated again to the flux quantum is too short in comparison to the interparticle separation -^ = A (which 
fitted to the B field orbits but not to 35 orbits). This is illustrated in the central panel of Fig. [2j Nevertheless, 
if tree-loop orbits are considered, then in flat geometry of 2D space, the external flux 3BA must be shared between 
three loops with the same surface A (i.e., BA for each loop). Thus, each loop accommodated to the flux quantum 
has the orbit with the surface A and gives the contribution BA to the flux, resulting in the total flux 3BA per 
particle, as needed—this is schematically illustrated in Fig. [2] (right). The size of A in the right panel is equal to A 
in the left panel, which means that the three-loop orbits fit to the interparticle separation equaled to A. 

As the braid group generators must be defined by the half of the cyclotron orbit (cf. Fig. [3]) thus the braid with 
one additional loop corresponds to the cyclotron orbits with three loops—such a generator has the form = (tTi)^ 
and the braid group generated by 6) , z = I,..., is the subgroup of the original full braid group. This subgroup 
is called the cyclotron braid subgroup and its IDURs define statistics of 2D charged particles at strong magnetic 
field corresponding to fractional filling z/ = | of the LLL giving rise to the explanation of the Laughlin statistics for 
FQHE at this fractional filling. The generalization to more loops attached to the braid generator one by one, results 
in double increase of loop number in multilooped cyclotron orbits and thereby in fractions i/ = p — odd integer. 
This approach successfully reproduces the hierarchy of the experimentally observed filling fraction corresponding to 
FQHE in the LLL and in the higher LLs for 2DEG in conventional semiconductors m- 

If ^ < A, i.e., when cyclotron orbits are larger that interparticle separation, (as in the right panel in Fig. [T|), some 
special commensurability circumstances important for braid group definition also occur. For the filling fractions when 
x-^ = A; X — integer, the cyclotron orbits fit to every cc-th particle separation (x-^ = j^), which also allows for the 
definition of the generators <Ti in the form of ordinary singlelooped braids, similarly as happens for the completely 
filled higher LLs. The related statistics is the same as for IQHE though at some fractional fillings of higher LLs as is 
demonstrated in Ref. in good cor resp ondence with experimental observations up to the three first LLs for 2DEG 
in conventional semiconductors [l,[20,l2y . This could happen only in such LL subbbands where the condition x-^ = A 
( X — integer) could be fulfilled. In the conventional semiconductor Hall systems it may happen only for n > 1 where 
n is the number enumerating LLs. Simultaneously, for zr > 0 too short cyclotron orbits can be encountered only 
close to subbands borders, because in the higher LLs the cyclotron singlelooped orbits are larger as accommodated to 
higher kinetic energy. This pushes FQHE (multiloop) features in higher LLs toward the edges of subbands in LLs with 
n > 0, whereas in central regions of these subbbands the new fractional features occur related with the singleooped 
cyclotron orbits and thus with IQHE-type of correlations but at fractional filling ratios. These correlated states are 
referred as to FQHE(singleloop). The quantization of the transverse resistance Rxy related to these fractional filling 
ratios of higher LLs z/ = ^ is as for FQHE ^, but the correlations of Laughlin type are with the exponent p = 1 
displaying singlelooped braid exchanges like in IQHE. The number of these new fractional filling ratios grows as 2n 
with the LL number. 
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III. COMMENSURABILITY INDUCED HIERARCHY OF FRACTIONAL FILLINGS FOR GRAPHENE 

In graphene one deals with the relativistic version of LLs [1, . This is due to specific band structure in this 

material being a gapless semiconductor with points K and K' on the border of the hexagonal Brillouin zone where the 
valence and conduction bands meet together Q . Thus the low energy particle-hole excitations can be described by the 
effective Dirac Hamiltonian corresponding to cone-shape of both bands close to the meeting points. The quantization 
due presence of the magnetic field has the form of degenerated LI^ though with spectrum not equidistant as for 
ordinary 2DEG but enumerated by ^/n, n is the number of the LL . This form of the selfenergies results from the 
linear in momentum Hamiltonian close to K and K' points, whereas the degeneration of each LL subband is the same 
as for 2DEG and equals to 7 ^^- The number of subbbands per each LL (each n) is here 4. This corresponds to the 
ordinary Zeeman spin-splitting and to the so called valley splitting (expressed often in terms of a pseudospin) due to 
doublet of inequivalent K points mixed with two sublattices for C atoms in crystal net of graphene [ 2 ^ . Taking into 
account that the Zeeman splitting in graphene is small @ and the valley splitting depending on the external magnetic 
field is small as well, the 4-fold approximate spin-valley additional degeneration is employed to determine the filling 
fractions for relativistic IQHE in graphene, in the following form v = 4(n-|- in pretty good correspondence with the 
experimental observations Q . The remarkable difference between this filling rate formula and that one for the ordinary 
2DEG (with factor 2 instead of 4, due to only approximate spin degeneration) is the presence of the overall shift by 
the factor 2 resulting from i in the above formula. This shift is due to the so-called Berry phase manifesting in the 
LLL for graphene and resulting in the sharing of states from the LLL between particles and holes from the conduction 
and valence band, respectively, at the zero energy level. Due to this feature the bottom of the LLL is shifted by 
2 (in terms of the filling factor) upward if counting only negatively charged carriers and oppositely—downward for 
positive holes. The corresponding filling rates for holes from the valence band are negative mirror refection of those 
for electrons from the conduction band. Note that changing between particles and holes can be easily achieved in 
graphene by the shift of the Fermi level around the Dirac point by application of a lateral voltage. 

The Berry phase contribution is referred sometimes as to the additional tt phase shift due to the chirality induced 
by the valley pseudospin when one adiabatically traverses with a selected particle a closed loop (e.g., along the 
semiclassical cyclotron loop) when the momentum made the 27r convolution. 

Similar analysis of the LLs can be done for the bilayer graphene [^ . . Due to off-diagonal interlayer hopping 

term the local Hamiltonian attains again the quadratic form with respect to the momentum. Thus the LL spectrum in 
bilayer graphene resembles that one for the ordinary 2DEG with four subbands for each level except for the LLL which 
is eight-fold degenerated. This extra degeneration of the LLL arose from the action of the square of an annihilation 
operator on oscillator states with n = 0 and n = 1 [22|, [ 2 ^ . Due to the sharing of the LLL between particles and 
holes, the bottom for uniformly charged carriers is located in the center of the 8 -fold quasi-degenerated LLL. This 
property is associated again with the Berry phase shift for chiral particles, though in the case of bilayer graphene 
it gives additional 27r phase shift. Therefore, steps of the relativistic IQHE are located in bilayer graphene again at 
intege^ filling rates of the subsequent LLs whereas for monolayer graphene were located at half-fillings of the LLs 


A. Monolayer graphene 

For magnetic field strong enough that v & (0,1) one deals with completely filled valence band hole states of the LLL 
and with the fractional filled subsequent particle states from the conduction band subband n = 0 , 2 f (2 indicates the 
valley pseudospin orientation and f indicates orientation of the ordinary spin). The degeneration of each subband is 
Nq = 7 ^^ and we have less than Nq of electrons, N < Nq for fractional fillings. 

The cyclotron orbits in the LLL must be accommodated to the bare kinetic energy T = hu}c{n+ i) where uic = 

The cyclotron orbits have the same size for all particles due to the flat band condition quenching the kinetic energy 
competition and resulting in the same averaged velocity and the same cyclotron orbit size for all particles at a 
perpendicular magnetic field presence (though quantumly, a velocity is not well defined since its coordinates do not 
commute). The cyclotron orbits restrict the topology of all trajectories uniformly in 2D, thus restrict the braid group 
structure despite of particularities of interaction and other fields including the crystal field. Therefore for graphene 
the cyclotron orbit structure is governed by ordinary (the same as for 2DEG) Landau levels restrictions even though 
a specific band structure with Dirac points being the result of the crystal field highly modifies LLs but not in terms 
of the bare kinetic energy. The particularities of graphene quantum dynamics induced by the electric-type interaction 
are included in the ordinary part of the Feynman path integral, whereas the additional summation over topologically 
nonequivalent trajectory classes concerns the braid group structure the same as for 2DEG upon the magnetic field. 
The difference between the conventional 2DEG systems and graphene will be related in this regard with distinct 
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number of LL subbands in graphene comparing to the conventional semiconductor 2DEG and in the Barry phase 
induced shift of fillings in the LLL. 

Thus, the cyclotron orbit size in the subband n = 0, 2 f equals to, = ■^, where S is the sample surface. As this 
orbit size is lower than the interparticle separation ^ (as N < Nq), the multilooped braid structure is necessary. From 
the commensurability condition one finds v = {q — odd integer to protect the braid structure). For 

holes in this subband one can expect the symmetric filling ratios v = 1 — ^ ■ Similarly as for the ordinary 2DEG one can 
generalize this simple series by assumption that the last loop of the multilooped cyclotron orbit can be commensurate 
with the interparticle separation as for some other filling ratio expressed by /, whereas the former loops take away an 
integer number of flux quanta. In this way one can obtain the hierarchy of fillings for FQHE in this subband of the 
LLL, V = j-i 1 v = 1 — ; where ^ = 1, 2,... or equals to other fractional filling factor and minus in the 

denominators indicates possibility of the eight-figure orientation of the last loop with respect to the antecedent one. 
The Hall metal states can be characterized by the limit 1 —>■ oo in the above formula (i.e., for zeroth flux taken away 
by the last loop as for ordinary fermions without the magnetic field, which was the case in the Hall metal archetype 
for ^ i in the conventional 2DEG), which gives the hierarchy for the Hall metal states, v = v = 1 — To 
account for the Berry phase anomaly in graphene the overall shift of v by —2 can be performed, but we use here the 
net filling fractions. 

For the completely filled subband n = 0,2 f, i.e., for v = 1, one arrives at IQHE. For lower magnetic field (or larger 
number of electrons), when three first subband of the LLL are filled and in the last subband of LLL, the cyclotron 
orbit size is still lower than the interparticle separation (because N — Nq < Nq). Thus multilooped structure 

is repeated from the previous subband. This results with the same FQHE hierarchy as for antecedent subband, only 
shifted ahead by 1. For the completely filled LLL (i.e, for completely filled all its subbands) one deals with IQHE 
according to its main-line = 4(n -I- ^). 

Similarly one can consider fillings of the second LL with n = 1. This level also has four subbands, but in this level the 
bare kinetic energy is equaled to and the related cyclotron orbit size ^ = ■^. For N G (2A^o, 3A^o] we deal with 
graduate filling of n = 1,1 f subband. Gyclotron orbits of size ^ are comparing here to interparticle separation scale 
N- 2 No • Only for small number of electrons in this subband one deals with the multilooped structure (corresponding 
to the inequality < jv-^Afo '''^^en q^ = jv-'ljVo ’ q~odd integer, which gives main series for FQHE (multiloop) in 
this subband, n = Similarly as before the complete related hierarchy reads, i/ = 2-1- n = 4— 

with the Hall metal hierarchy in the limit I —?► oo. These series are located closer to the subband edges, whereas in 
the center of this subband the other commensurability conditions are possible. When ^ ~ 2, 3 we 

get n = ^, |,3, respectively, corresponding to singlelooped cyclotron orbits similar as for IQHE. Thus, for i/ = |, | 
one deals with FQHE(singleloop)—the new Hall feature manifesting itself only in higher LLs, where cyclotron orbits 
may be larger than the interparticle separation. At the special case arrive at iz = | with the 

paired particles (pairing does not change the cyclotron radius but twice diminishes the carrier number to ). 

The following subband are filled with electrons upon the similar scheme. For the subband n = 1,1 i the cyclotron 
size is whereas the interparticle distances are measured with the plaque where N G (3iVo,4iVo). The 

commensurability condition results in the main series for FQHE (multiloop), n = 3 + ■^, which can be 

developed to the full hierarchy similarly as described above. The condition with a; = 1,2,3 results in 

fractions with singleloop correlations of FQHE(singleloop)-type for n = ^ and IQHE for iz = 4, correspondingly, 

whereas a paired state can be realized at zz = |. 

B. Bilayer graphene 

The special topology of the bilayer 2DEG structure creates opportunity to verify the braid group based concept of 
the commensurability of cyclotron orbits size with the interparticle spacing in planar system of interacting particles. 
The interaction preserves particles do not approach one to another one closer than particle separation resulting from 
the planar density. Only when the cyclotron orbits fit accurately to this spacing the mutual classical exchange of 
neighboring particles is possible at presence of the perpendicular magnetic field. This is caused by the fact that 
the trajectories of particles are governed by the cyclotron orbital action of the magnetic field and despite of the 
interaction particularities the geometry of free particle cyclotron orbits allow for topological classification of possible 
or not possible trajectories creating the fundamental tti group of the configuration space for the system of N particles 
on the plane, i.e., the full braid group for the plane |13j . In 2D all cyclotron orbits are planar and in the case of 
uniformly distributed particles with the same velocities exposed to the perpendicular magnetic field they either admit 
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topologically exchanges of neighboring particles or not depending on the commensurability of the cyclotron orbits 
with the interparticle separation. If the planar size of the cyclotron orbits A equals to the plane fraction per single 
particle A = ^ (where S is the surface field of the sample and N is the number of particles), the mutual exchanges of 
particles are possible and the ordinary full braid group 7ri(<I>) can be defined. Otherwise, when ^ the exchanges 
of particles are impossible along ordinary cyclotron trajectories, as they are too short to merge neighbors on the 
plane. When ^ ^ the interchanges of 2D particles are also impossible along cyclotron trajectories, because such 

an exchange do not allow for conservation of uniform particle distribution with the constant interparticle spacing. 
These three distinct commensurability situations were depicted in Fig. [T] Topological conditioning of interparticle 
exchanges in the uniformly distributed planar systems of charged particles plays the fundamental role in discrimination 
of trajectories from the full braid group which are not admitted at the presence of the magnetic field as do not merge 
particles. Miss mash between size of cyclotron orbits and interparticle spacing precludes the definition of the full 
braid group generators di corresponding to the interchange of i-th and i + 1-th particles (upon a selected certain 
enumeration of particles, arbitrary in general due to particle indistinguishability). When the braid group cannot be 
defined the collective state cannot be created as the statistics of particles cannot be determined. In order to establish 
the statistics of quantum particles, the related braid group must be defined and the statistics is governed by the 
selected IDUR of this group [T^ . 

For the left panel case in Fig. [T] the full braid group is well defined and all statistics typical for 2D are available, 
i.e., according to the form of IDURs for the 7ri(<I>), M = R^, e®“ with a € [0,27r) corresponding to fermions for 
a = TT, to bosons for a = 0 and to anyons for other a. If, however, -^ > A (as in the central panel in Fig. [1}, the 
ordinary generators ai cannot be defined. If one removes from the full braid group all impossible singlelooped-related 
trajectories corresponding to generators of this group, which cannot be defined as the cyclotron orbits are too small to 
match particles, one can observe that the remaining braids with additional loops can match particles on 2D separated 
larger than the singlelooped cyclotron orbit reaches. These braids are related to multilooped cyclotron orbits with 
p loops, i.e., when pA = with p odd integer 0- This property follows from the fact that multilooped cyclotron 
orbits on the plane must be larger than the singlelooped orbits for the same magnetic field. The external magnetic 
field flux passing through the multilooped 2D cyclotron orbit is the same as that one passing through the singlelooped 
orbit and in the former case, per each loop falls only fraction of the total flux. It means that the size of each loop 
grows to fit to the fraction of the flux as illustrated in Fig. [2] Repeating again the related argumentation we show the 
cyclotron orbit at magnetic field B in Fig. [2] (left) as accommodated to the quantum of the magnetic field flux, i.e., 
BA = because the cyclotron orbit A fits to S —the sample area, N —the number of particles, in the case of 
the completely filled LLL. At p-times larger field, pB, the cyclotron orbit accommodated again to the flux quantum 
is too short in comparison to the interparticle separation as is illustrated in the central panel of Fig. [^for p = 3. 
If p-looped orbits are considered, then in of 2D space, the external flux pBA must be shared between p loops with 
the same surface A (i.e., BA for each loop). Thus, each loop accommodated to the flux quantum ^ has the orbit 
with the surface A, resulting in the total flux pBA per particle, as is illustrated in Fig. [2] (right). The size of A in 
the right panel is equal to A in the left panel, which means that the p-looped orbits fit to the interparticle separation 
defined by A though the singlelooped not. This property attributed exclusively to the exact 2D topology allows for 
explanation of FQHE and related exotic Laughlin correlations [l^ . 



FIG. 1: Schematic demonstration that commensurability (left) of cyclotron orbit with interparticle separation satisfies topology 
requirements for braid interchanges in equidistantly uniformly distributed 2D particles; for smaller cyclotron radii particles 
cannot be matched (center), for larger ones the interparticle distance cannot be conserved (right) 

The bilayer graphene is, however, not strictly two dimensional and for the bilayer graphene the topological situation 
changes considerably. Two sheets of the graphene plane lie in close distance with the hopping constant allowing 
changing of electron positions between the planes. Thus, we deal here with double number of electrons residing on 
two-sheet structure instead of the single sheet as was the case for the monolayer graphene. 

All above described requirements to fulfill commensurability condition in order to define the related braid group 
describing the correlated multiparticle state are in charge also for the bilayer graphene, with a single distinction with 
respect to the monolayer case. The doublelooped cyclotron orbits may have in bilayer graphene the same size as the 
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A B = hc/e 



A/3 3B = hc/e 



A 3B = 3 hc/e 



FIG. 2: Schematic illustration of cyclotron orbit enhancement in 2D due to multi-loop trajectory structure (third dimension 
added for visual clarity) 


single exchange - braids 


double exchange - cyclotron orbits 



FIG. 3: The braid generator at corresponds to the single exchange of particles (left), the cyclotron orbit (relative) corresponds 
to the double exchange (right), for v = 1 when singlelooped cyclotron trajectory reaches neighboring particles, Rc = Ro (a) 
braid generator erf for v ~ \ with additional two loops needed for Rc = Ro (i?c-cyclotron radius, 2iio-particle separation, i.e., 
-nRl = itRI = I) (b) 


singlelooped orbit. This follows from the fact that the second loop can be located in the opposite sheet of graphene 
than the first one and the external field passing through such doublelooped orbit is twice larger than the flux passing 
through the singlelooped orbit. Each loop has in this case the separate own surface in contrary to the multilooped 
cyclotron orbit located on the purely 2D plane which forced each loop to take away only some fraction of the total 
flux because all these loops share the same surface in 2D. Taking into account that in the bilayer system loops of the 
multilooped orbit may be located partly in both 2D sheets, the contribution of the one loop must be avoided whereas 
the remaining loops must share the same flux as passing through a singlelooped orbit, independently how loops are 
distributed between two sheets. Thus, one can write out the commensurability condition in the bilayer graphene for 
the case of too short singlelooped cyclotron orbits in the following form (for concreteness in the subband n = 0, 2 f of 
the LLL—the first particle-type subband of the LLL), 


he 

eB 


= — < 
Ar_ ^ 


N-Nq^ 

(p-l)§ = (p-l)f = 
= JL = 1 _iii 


No 


p-1 


2’ 4’ 6’ ■ 


N-Nq ’ 


(3) 


where, N is the total number of particles in both graphene sheets, Nq is the degeneration counted for both sheets 
together, S is the surface of the sample (the surface of the single sheet) p is an odd integer to assure that the half of 
the cyclotron orbit defines the braid. 

The factor p — 1 in l.h.s. of Eq. ([3]) is caused by the fact that in enlarging of the effective cyclotron orbits participate 
only orbits from the ideally 2D sheet of bilayer graphene (no matter in which are located doubling loops) with the 
exception of a single orbit which is located in opposite sheet to the first one. This sole loop contributes to the total 
flux with the additional flux quantum due to its own surface and this loop must be omitted. The next orbits must 
duplicate the former ones (in fact two) without rising the surface and no matter in which sheet are they located, 
because in both they will duplicate loops already present there. Thus, in the enhancement of the effective p-looped 
















cyclotron orbit take part only p — 1 loops. 



paired 

n=1 


FQHE(singleloop) 



net filling factor = N/N„ - 2 
paired 


FQHE(singleioop) 



net filling factor = N/N„ - 6 

FIG. 4: Evolution of fractional filling hierarchy in three first LLs of monolayer graphene; for each LL the first particle subband 
is illustrated, the next subbands in each LL repeat the hierarchy from the first one. Different types of ordering are indicated 
with spikes of various height. Series for FQHE (ordinary—multiloop), FQHE(singleloop), Hall metal and paired state are 
displayed acc. to the hierarchy described in Tab. |T]with q = 3 — 9, 1 = 1 — 20; only a few selected ratios from these series are 
explicitly written out. 

Let us emphasize that for such multilooped orbits the total number of loops still is p —thus the generators of the 
corresponding cyclotron subgroup are of the form = af, resulting in the Laughlin correlations with the p exponent 
for the Jastrow polynomial. But due to the distinct commensurability of orbits with interparticle separation ([3]) the 
related filling fraction are iz = in the first particle-type subband of the LLL, i.e., the subband n = 0, l,t. This 
even denominator main series of the FQHE hierarchy for bilayer graphene pretty well coincides with the experimental 
observations [H. 

For holes in this subband (let us emphasize that these holes are not holes from the valence band but correspond to 
unfilled states in the almost filled subband of particle type) one can write, iz = 1 — whereas the generalization to 

the full hierarchy of FQHE in this subband attains the form, i' = ;^p_ 2 )±i ’ = ^ — ;(p_ 2 )±i ’ ^ corresponds to 
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some filling factor for other correlated Hall state, including completely filled LLs with IQHE. In the next subbands 
of the LLL, n = 0, 2 | (assuming tht this subband succeeds the former one), the hierarchy is identical only shifted by 
one, because commensurability condition has the same form for all subbands with the same n due to the same size of 
the cyclotron orbits. 

Some novelty occurs, however, in the following two subbands of the LLL, n = 1,2 f and n = 1,2 The FQHE 
main series in the first of these subbands of the LLL, n = 1, 2 f, has the form, 


3/ic _ 35 ^ S 

eB ~ No ^ N-2 Nq ’ 

{P — 1 ) 3 ^ = (j> — 1 )|^ = N-2No ’ 

=2 + ^ = 2+i2+i,4- 


18’ 


(4) 


The generalization of this main series for holes in the subband and to the full FQHE hierarchy in this subband is as 
follows: for subband holes (i.e., for not completely filled subband), z/ = 3 — and for the full FQHE hierarchy in 

this subband, z/ = 2 + 3 ;(pJ 2 )±i ’ ^ = 3 — 3 ;(pJ 2 )±i Hall metal hierarchy in the limit I ^ oo). 


LL subb. 

FQHE(singleloop), paired, IQHE 

FQHE(multiloop) {q — odd, Z = 1, 2,3,...) 

Hall metal 

n = 0, 2 t 

V=1 

j/ — 1 I * i * 

^ = “S'! 1-S' 

q-l’ 9-1 

n= l,lt 

V — (| paired), 2,3 

,, _0|1q ‘ 3 ' 

^ 3Z(<?-1)±1’'^ 3/(q-1)±1 

2+3(^^^,3 

n = 2,1 t 

ai aa a4 / la nnirpd\ fi 7 

fi 1 (S 1 * 7 7 * 

?■ fi 1 7 


^ ^ 1 ' 51(q-1)±1^‘ 5q’' 5l(q-l)±l 

° 1 5(9-1)’ ‘ 5(9-1) 


TABLE I: LL filling factors for FQHE determined by commensurability condition {paired indicates condensate of electron 
pairs), for the first particle subband in each of the three first LLs (n = 0,1, 2) for monolayer graphene 


LL subb. 

FQHE(singleloop), paired, 

IQHE, 1/ = 

FQHE(multiloop) {q — odd, Z = 1, 2, 3,...), z^ = 

Hall metal, v = 

n = 0,2t (0) 

0,1 

1 1 i i 1 1 

(9-1) ’J- (9-1), Z(9-2)±1> J- Z(9-2)±1 

-J— 1_ L- 

q-2 5 q-2 

n = l,2t(0) 

5 . (| paired), 2,3 

0|4q loi !■ O 1- 

3fq-l)’“ ' 3l(q-2)±l^^ 3l(q-2)±l 

3(9-2)’'^ 3(9-2) 

n = 2,lt (1) 

T’ ¥’ ¥’ T’ (I paired), 4,5 

T+ 4 4- ‘ 5 1 r 1 

' 5l(q-2)±l^^ 5(q-l)^^ 5l(q-2)±l 

a 1 1 1 

^+5(9-2)’° 5(9-2) 


TABLE II: LL filling factors for FQHE determined by commensurability condition {paired indicates condensate of electron 
pairs), for the first particle subband in each of the two first LLs (n = 0,1 for the extra degenerated LLL and n = 2 for the first 
LL above the LLL) for bilayer graphene 


LL subb. 

FQHE(smgleloop), paired, 

IQHE, z/ = 

FQHE(multiloop) (g — odd, Z = 1, 2, 3, ...), z/ = 

Hall metal, v = 

n = 0,2t 

1 

1 1 i i 1 1 

(9-1) ’J- (,1), i(„2)±l> J- Z{9-2)±1 

-J— 1 — ZEI 

q-2 ’ ^ q-2 

n = l,2t 

|, §, (| paired), 1,2 

1 1 1 O ^11 ‘ O Z 

3fQ-l)’-^ ' 3Uq-2)±l’" 3l(q-2)±l 

^+3(9-2)’“ 3(9-2) 

n = l,2t 

(I paired), 1 

1 / 1 1 1 

3(q-l) ’ 3Z(q-2)±1 ’ 3(q-l) ’ ^ 3l(q-2)±l 

3(9-2)3(9-2) 

n = 0,2t 

1,2 

1+^1+ ‘ 2 ^2 ‘ 

Z(q-l)±l’" q-2’"- Z(q-2)±1 

1 j_ L- 2 _ — 

^ q-2 ’ ^ q-2 


TABLE HI: Comparison of filling hierarchy in the LLL level in bilayer graphene for two mutually inverted successions of two 
lowest subbands: rz = 0, 2 b, u = 1, 2, t (upper) and n = 1, 2 b, u = 0, 2 b (lower) 

Nevertheless, in the subband n = 1, 2 b of the LLL some new commensurability opportunity occurs: 
for X = 1,2,3, which gives fillings ratios v = |,|,3, correspondingly. All this ratios are related with singlelooped 
cyclotron trajectories, thus with singleloop correlations similar as for IQHE (though for not integer filling ratios). 
This new Hall feature, typical for LLs with n > 1, we called as FQHE (singleloop). Moreover, for x = 1.5 one can 
consider twice diminishing of particle number {N — 2Nq)/2 due to the pairing, which gives perfect commensurability 
of cyclotron orbits of pairs with the separation of the particle pairs at z/ = |. 

The last subband n = 1, 2 b in the LLL in bilayer graphene is filled with electrons in the similar manner because for 
both subbands with n = 1 the cyclotron orbits have the same size. Thus the hierarchy of fractional filling for the last 
subband in the LLL is shifted by 1 from the antecedent subband without any modification. The situation changes, 
however, in the next LL (the first one above the LLL). In the hrst such LL (with n = 2) the cyclotron orbits suited 
to commensurability condition are determined by the bare kinetic energy for rz = 2, and the corresponding cyclotron 
orbit size is equaled to, The similar analysis as in the previous LL, gives here the main series and the full 
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hierarchy for FQHE(multiloop) in the subband n = 2, If, v = A + u = A + 5 ;(pJ 2 )±i ’ respectively (inclusion 

of subband holes resolves itself to the substitution of 4+ by 5— in both above formulae). Similarly as previously, the 
limit I —>■ oo gives the Hall metal hierarchy. The difference in comparison to the previous LL consists here also in 
the presence of four (instead two) satellite FQHE (singleloop) states symmetrically located around the central paired 
state. In the subband n = 2, l,t the satellite states occur at ^ and the central paired state at i/ = |. 

This hierarchy is repeated in all 4 subbands of the first LL. 

The evolution of the fractional filling hierarchy of subsequent LLs is illustrated in FiglU and Fig. [5l for monolayer 
and bilayer graphene, correspondingly and is summarized in Tab. Uand Tab. Ih] 

For bilayer graphene the degeneration of n = 0 and n = 1 states results in 8-fold degeneration of the LLL doubling 4- 
fold spin-valley degeneration. The degeneration is not exact and with rising magnetic field amplitude both the Zeeman 
splitting and the valley splitting grows. Stress, deformation and structure imperfections also cause the increase of 
the valley splitting. Inclusion of the interaction plays a similar role. Coulomb interaction causes mixing of n = 0,1 
states lifting their degeneration. Especially interesting is such a degeneration lifting which admits inverted order of 
fillings of LLL subbands with distinct n. The inversion of orders n = 0,1 to n = 1,0 affects the filling ratio hierarchy. 
Assuming that the LLL subbands with n = 1 is earlier filled than the n = 0 subband, we get the following hierarchy 
for the first subband n = 1, 2 f: multilooped orbits for v = 2 i(p-i)±i ^ = 1 — 3 ;(p-i)±i ; siglelooped orbits for = i, | 

and paired state for = ^- Assuming the next subband, n = 0, 2, f, we get the hierarchy of fillings for this subband 
in the form: multilooped orbits for = 1+ ;(p_i) j-i > = 2 — and no singlelooped orbits. The comparison of 

reverted orderings of two first LLL subbands is summarized in Tab. IIIII 

One can consider also the situation in the LLL of bilayer graphene, when the degeneration of n = 0,1 states is 
lifted is such a way that both levels cross at certain filling factor zz* < 1 (cf. [IJ where mixing between n = 0,1 
states is numerically analyzed for small models on torus or sphere). Let us assume for a model, that first is filled the 
n = 1 subband (n = 1, 2 f) up to zz*. At this filling the subband n = 1, 2 f crosses with the subband n = 0, 2 f and 
the latter is filled for 1 -|- zz* > v > v*. The related hierarchy of fractional filings looks like an ordinary filling of the 
subband n = 1, 2 f, however, with an insertion of n = 0, 2 f. Depending on the value of zz* the various patterns are 
achievable by simple combination of hierarchy patterns listed in Tab. IIIII (including also inverted ordering of n = 1 
and n = 0 subbands). 


IV. COMPARISON WITH EXPERIMENT 

Searching states related to FQHE in the case of Hall measurements in graphene is particularly challenging since 
the different ’relativistic’ structure of LLs more complicated in comparison to the conventional semiconductor 2DEG 
Hall physics. Moreover the filling factor can be changed in graphene both by the external magnetic field and by 
the particle concentration via shifting of the Fermi level near the Dirac point by application of the lateral voltage. 
Due to spin-valley degeneration and the Berry phase contribution related to chiral valley pseudospin, the IQHE is 
observed in graphene for fillings zz = 4(n-|- i) = 2, 6,10,14,.. for particles from the conduction band and for the same 
but negative factors for holes from the valence band. Despite using very strong m^netic fields (up to 45 T), FQHE 
was, however, not detected in graphene samples deposited on a substrate of Si02 [7|. Instead, at so strong magnetic 
fields the emergence of additional plateaus of IQHE has been observed for the fillings zz = 0,±1,±4, indicating the 
elimination of spin-valley degeneration, as a result of increasing mass of Dirac fermions Q- Only after mastering 
the novel technology of the so-called suspended ultrasmall graphene scrapings with extreme purity and high mobility 
of carriers above 200000 cm^V“^s“^ (high mobility is necessary to observe FQHE also in the case of conventional 
semiconductor 2D hetero-structures, which may be related to multi-looped quasi-classical cyclotron movement of wave 
packets in the case of multi-looped braids associated with FQHE 1^: note that in semiconductor 2D heterostructures 
carrier mobility reaches even higher values of millions cm^V“^s“^ |26l|L it was possible to observe FQHE in graphene 
at net fillings zz = 1/3 and —1/3 (the latter for holes, at the opposite polarization of the gat e voltage, which determines 
position of the Fermi level, either in the conduction band, or in the valence band) [^l28|. Both these papers report 
the observation of FQHE in graphene for medium strong magnetic fields. The paper [^, in a field of 14 T, for 
electron concentration of 10^^/cm^ and the paper [2^, in a field of 2 T, but for a concentration level smaller by one 
order of magnitude (10^°cm“^ and the mobility of 200000 cm^V“^s“^) as shown in Fig. |6l 

FQHE in suspended graphene is observed at the temperatures around 10 K [2^, and even higher (up to 20 K) [30l| . 
which seems to be related with the stronger electric interaction in view of lack, in the case of suspended samples, of 
a dielectric substrate (with the dielectric constant in case of Si02, ^ 3.9) and, on the other hand, with very high 
cyclotron energy in g:rap hene (i.e., large energy gap between incompressible states). 

In the papers l32| it has also been demonstrated the competition between the FQHE state with the insulator 
state near the Dirac point, corresponding to a rapidly decreasing concentration—Fig. [T] 
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FIG. 5: Evolution of fractional filling hierarchy in two first LLs of bilayer graphene; for the LLL two subbands with n = 0 
and n = 1 are illustrated. Different types of ordering are indicated with spikes of various height. Series for FQHE (ordinary— 
multiloop), FQHE(singleloop), Hall metal and paired state are displayed acc. to the hierarchy described in Tab. [U with 
g = 3 — 9, / = 1 — 10; only a few selected ratios from these series are explicitly written out. 


From the perspective of cyclotron groups, experimental results on FQHE in graphene p7l - [^ Is^ seem to be 
compliant with the expectations of the braid description. In the case of graphene, the specific band structure with 
conical Dirac bands leads to simultaneous participation (in Dirac point) of both bands—of holes and of electrons, 
which combined with the massless character of Dirac fermions manifests itself through an anomalous ’relativistic’ IQHE 
Sim, HI. Controlling lateral gate voltage (within the range ca. 5 — 60 V [Il,i[23) allows regulating the density 
of carriers at a constant magnetic field. One should therefore expect that at relatively small densities of carriers 
(electrons, or symmetrical holes at reverse voltage polarization), the cyclotron orbits will be too short to prevent 
braid exchanges of particles at a sufficiently strong magnetic field—although weaker for smaller concentrations—and 
experimental observations exactly support that j27l l28l|. Eor low concentration, while closing on the Dirac point, one 
may expect that too strong fields would exceed the stability threshold of the EQHE state in competition with the 
Wigner crystal (assuming a similar character of this competition in the case of massless Dirac fermions in reference 
to conventional semiconductor 2D structures) and that corresponds to the emergence of the insulating state near the 
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H {10'^ cm ^g{V) 

FIG. 6: a) FQHE observation in suspended graphene for the filling 0.3 (1/3) in a field of 12-14 T with the concentration of 
lO^^cm”^ and the mobility of 250000 cm^V~^s~^, b) FQHE singularities in suspended graphene for the filling | in a field of 
2-12 T with the concentration of 10^°cm“^ and the mobility of 200000 cm^V“^s“^ (after [^. 




FIG. 7: a) The emergence of an insulator state accompanying the increase in the strength of a magnetic field around the 
Dirac point, b) competition between FQHE and the insulator state for the filling —1/3: annealing removes pollution—enhances 
mobility and provides conditions for the emergence of plateau for FQHE (after [2^) 


Dirac point in a sufficiently large magnetic field [s^ . In the case of the hexagonal structure of graphene, electron (or 
hole) Wigner crystallization [36l | may exhibit interference between the triangular crystal sublattices, and including of 
the resonance (hopping) between these two sublattices may cause blurring the sharp transition to the insulator state, 
which seems compliant with observations (Fig. [7]). 

The recent progress in the experiment allows also for observation of FQHE in graphene on the crystal substrate of 
boron nitride {BN) in large magnetic fields of order of 40 T (remarkably, FQHE features were noticed in this case up 
to z/ = 4) 0. 

The mobility of carriers in graphene is lower than in traditional 2DEG, but taking into account that the carrier 
concentration in graphene can be lower in comparison to semiconductor heterostructure , the corresponding 

mean free path in both cases well exceeds the sample dimension (of pLva. order, as the mobility is proportional to the 
concentration and to the mean free path of carriers). 

Energy gaps protecting incompressible FQHE states are in graphene larger than in traditional semiconductor 
materials, reaching an order of 16 K (at z/ = | and B = 35 T), which is referred to Dirac massless character of 
carriers. In conventional semiconductor heterostractures the corresponding gaps are much lower and the observed 
FQHE stability with temperature is much more fragile. 

The recent development in experiments with monolayer graphene on BN substrate ai and with suspended small 
sheets ii allowed for observation of more and more Hall features at fractional fillings of subsequent subbands of two 
first LLs. While the sequence of fillings in the lowest subband of the LLL hts well to CF predictions (including CFs 
with two and four flux quanta attached), an explanation of the filling structure of next subbands strongly deviates 
from this simple picture. The pattern of filling rates repeated in the subbands of the first LL has apparently nothing 
in common with the CF concept Sii- This phenomenon is referred in these papers to the various scenarios of 
breaking of the approximate SU(4) spin-valley symmetry in graphene. Because of the smallness of the Zeeman splitting 
Ez in comparison to the Coulomb energy in graphene, Ez/Ec ^ O.Ole, similarly as of lattice scale in comparison 

to the magnetic length, a/ls ^ 0.06 (where = y and e is the dielectric susceptibility ~ 3.2) Q, the subbands 
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FIG. 8: Fractional quantum Hall effect for graphene on BN, a,b. Magnetoresistance (left axis) and Hall conductivity (right 
axis) in the n = 0 and n = 1 Landau levels at B = 35 T and temperature ~ 0.3 K (after i)- All filling ratios indicated in blue 
agree with the hierarchy given in Tab. [H 



Vg (Volts) 

FIG. 9: Magnetoresistance (left axis) and Hall resistance for graphene on BN (right axis) versus gate voltage acquired at B = 
35 T. Inset shows SdH oscillations at Vg = -18.5 V (after i). All filling ratios indicated in the figure (in blue) agree with the 
hierarchy as given in Tab. |T] 


which differ with spin and vaiiey-pseudospin orientation are cioseiy iocated and can be regarded as approximateiy 
degenerated. This SU(4) spin-vaiiey symmetry can be next broken by various factors and one can search arguments 
for unusuai fiiiing ratios hierarchy in reiated symmetry breaking and phase-iike transitions. Despite of many reiated 
ideas no fuiiy consistent picture is attained in this way as of yet. 

If one compares the experimentaiiy observed fractions for characteristic FQHE features in iongitudinai and Haii 
conductivities measured on variety of sampies of graphene with the pattern of fractionai hierarchy for two iowest LL, 
as iiiustrated in Fig. 01 one notices the coincidence of this hierarchy with the measured data. All fractions found 
experimentally can be reproduced by this hierarchy (cf. Tab. IJ). From this comparison it is visible why the CFs are 
efficient only in the LLL. This is linked with the fact that exclusively in the LLL cyclotron orbits are always shorter 
than the interparticle separation and additional loops are necessary. These loops can be modeled by fictitious held 
hux quanta attached to CFs. Though the analogy to additional loops is not exact, it allows to get at least the similar 
main line of the filling hierarchy in the LLL as that one given by the cominensurability condition. The usefulness of 
the CF model dramatically breaks down, however, in higher LLs because starting from the hrst LL the multilooped 
cominensurability is needed only close to the subband borders, whereas the central regions of all subbands of the 
hrst LL are occupied by doublets of hlling factors, (|, |), (^, ^), (^, ^), (^, ^), corresponding to singlelooped 
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commensurability condition, not allowing for CF modeling, but visible in experiments as FQHE(singleloop) [H|. The 
number of centrally located filling rates for FQHE(singleloop) grows next with the LL number as 2n. The repeating 
doublet of filling ratios for n = 1 is noticeably e.g., in Fig. |9]and in more accurate measurements in suspended samples 
ii besides of those on the BN substrate Q. 

The hierarchy induced by the commensurability condition reproduces the position of other observed features in two 
lowest LLs of the monolayer graphene. The elongated plateaus at border of subbands with IQHE seem to overwhelm 
also minima related by closely to border located FQHE rates being in this way out of the experimental resolution. 



FIG. 10: Observation of FQHE at T = 0.25 K in bilayer suspended graphene, a) magneto-resistance Rxx (blue curve) and R^y 
(black curve) at the lateral voltage —27 V (the inset shows an image of the actual device covered with a resist mask). Zoom-in 
on the data shown in a) to better illustrate the features associated with the —1/2 and —4/3 FQHE states, respectively; the 
features in the range marked by the yellow diamond are affected by remnant disorder (after i)- In red is added fitting with 
the hierarchy given in Tab. El 

The most convincing evidence supporting the correctness of the commensurability condition is, however, the coinci¬ 
dence of the related predictions with the experimental observations in the bilayer graphene. The pronounced feature 
of these observations is the occurrence of primary fractional features with even denominators also in the lowest sub¬ 
band of the LLL in bilayer graphene, oppositely to monolayer one. [I|. The commensurability condition for bilayer 
graphene reproduces perfectly the observed experimentally hierarchy—cf. Fig. 1101 as is illustrated in Fig. [S] and 
summarized in Tab. m 

Note finally that the FQHE hierarchy in bilayer systems with the characteristic even denominators holds also for 
bilayer 2DEG conventional Hall setups and indeed the = 5 state has been discovered there [H, . 


V. CONCLUSION 

The condition for commensurability of cyclotron orbits building cyclotron braid subgroups with interparticle spacing 
in uniform 2D charged systems is formulated in order to verify possibility of arrangement of correlated multiparticle 
states in planar quantum Hall systems. Using this commensurability condition based on the braid group approach to 
statistics of many particle systems the hierarchy of fractional fillings for LLs in graphene was determined. The FQHE 
evolution with growing number of LL has been described. The new opportunities for commensurability in higher 
LLs were established leading to different that ordinary FQHE(multiloop) correlated states at fractional fillings of LLs 
starting from the first one have been identified (they are referred to FQHE(singleloop) as the related correlations 
are described by singlelooped braids). Both the monolayer and bilayer graphene were considered and the essential 
difference of related hierarchy structures has been demonstrated and described. The even denominator main line of 
the fractional filling hierarchy in bilayer graphene is found in agreement with the experimental observations. The 
presented hierarchy for the monolayer and bilayer graphene found the confirmation in available experimental data for 
graphene on BN substrate as well as for suspended samples including bilayer suspended sheets for fillings of related 
spin-valley subbands from the zeroth and the first Landau levels. 
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